Cavendish-HEP-03/22 
LBNL-53701 



Top quark production at the Tevatron at NNLO 1 



Nikolaos Kidonakis a and Ramona Vogt 6 

^ ' a Cavendish Laboratory, University of Cambridge, 

O ■ Madingley Road, Cambridge CBS OHE, UK 

O ; 

. b Nuclear Science Division, 

Qh! Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 

C/3 , an d 

Physics Department, 

University of California at Davis, Davis, CA 95616, USA 



Abstract 



> 

in 

. We present results for top quark production at the Tevatron including next-to-next- 

\ to-leading order (NNLO) soft-gluon corrections. We show the stability of the cross section 

with respect to kinematics choice and scale when the NNLO corrections are taken into 
O ■ account. 
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1 Introduction 



Improved theoretical calculations of the top quark production cross sections and differential 
distributions are becoming more important as the top is observed at Tevatron Run II, where 
increasingly accurate measurements of the top mass and cross section are expected. The current 
state-of-the-art in the theory of top hadroproduction is calculations of next-to-next-to-leading 
order (NNLO) soft-gluon corrections to the double differential cross section PUI2] from threshold 
resummation techniques [HHUE]. Near threshold there is limited phase space for the emission 
of real gluons so that soft-gluon corrections dominate the cross section. 

These soft corrections appear in the form of logarithmic "plus" distributions [ln z (x t h)/^th] + , 
with I < 2n — 1 for the ra-th order corrections in a s , where x t h is a kinematical variable 
that measures distance from threshold. Calculations for top quark production until recently 
had been done through next-to-next-to-leading logarithmic (NNLL) accuracy, i.e. including 
logarithms with Z = 3,2, 1 at NNLO [HllZj- This NNLO-NNLL calculation greatly diminished 
the factorization and renormalization scale dependence of the cross section. However, the 
dependence of the corrections on the choice of kinematics, single-particle-inclusive (1PI) or 
pair- invariant-mass (PIM), was substantial even near threshold. Away from threshold, where 
hard real gluonic radiation becomes non-negligible, the discrepancy between 1PI and PIM 
results is not surprising. However, near threshold the results should be the same if all the 
NNLO soft corrections are included. Thus subleading (beyond NNLL) contributions can still 
have an impact on the cross section. 

We have recently calculated next-to-next-to-next-to-leading logarithmic (NNNLL) terms at 
NNLO in Ref. [Tj, i.e. terms with I = 0, following the unified approach of Ref. [2]. We find 
that the inclusion of NNNLL terms indeed eliminates the kinematics ambiguity near threshold. 



2 Analytical results 

We first discuss the analytical expressions for the corrections. We begin with the next-to-leading 
order (NLO) corrections. In the MS scheme, the NLO soft and virtual corrections for qq —* it 
in 1PI kinematics are 
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where F5, 1PI is the Born term and the renormalization scale. Also c|^~ = ACp and 
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1 — Am 2 /s. For use below, we 

write = 7^5 - 2C F ln( M |/ S ). Finally, cj" = a^f +V 1PI /[(«,/vr)F| ^] where a% S+V 1PI 
denotes the 5(s 4 ) terms in Eq. (4.7) of Ref. 8\ with the definitions of t\ and u\ interchanged 
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with respect to that reference. We also define c\ F ^ = T^ + C F [— 3/2 + ln(t 1 -u 1 /m 4 )] ln(^i F /s) + 



{Po/2) ln(/i^/s), where T^I has no scale dependence. 

We further define the constants C2 = tt 2 /6, and Cs = 1.2020569 • • •, /3 = (IIC A - 2n/)/3, 
and if = CU(67/18 — 7r 2 /6) — 5n//9 where n/ is the number of light quark flavors. 

Following Ref. |2j we write the NNLO soft-plus-virtual corrections in IPI kinematics as 
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Here the term = C F C A (7( 3 /2 + 22( 2 /3 - 299/27) + n f C F (-4:( 2 /3 + 50/27) denotes a set 
of two- loop contributions universal for processes with qq initial states [2] . Process-dependent 

(2) 

two-loop corrections are not included in Q but their contribution is expected to be negligible. 
The virtual contribution R^ 1 is not fully known. However, we can determine certain terms in 
Rg? exactly. These exact terms involve the factorization and renormalization scales as well as 
the those terms (( terms) that arise from the inversion from moment to momentum space PQ. 

The analytical form of the NLO and NNLO corrections for the qq channel in PIM kinematics 
is similar, see Ref. [Tj. The expressions for the gg channel in both IPI and PIM kinematics are 
more involved since the gg color structure is more complex pQ. 



3 Numerical results 

We now study the numerical impact of the corrections. The total partonic cross section may be 
expressed in terms of dimensionless scaling functions f^' 1 ^ that depend only on 77 = s/Am 2 — 1 
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Here we have set /i = fi F = u R . 

In Fig. ^ we plot the fij' ^ scaling functions, the most important contributions at NNLO 
and independent of \x. To demonstrate the effect of adding successive subleading contributions, 
we show the NNLL results in the upper plots, the scaling functions through NNNLL in the 
middle plots, and the results with the NNNLL and virtual ( terms in the lower plots. 

We first discuss the results for the qq channel. We note that to NNLL, the two kinematics 
choices give rather different results, even at low rj. When the NNNLL terms are added the two 
results coincide near threshold. Adding the virtual ( terms resulting from inversion improves 
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Figure 1: The ' scaling functions in the MS scheme for the qq (left) and gg (right) channels. 
The top plots show the NNLL result, the center plots the NNNLL results and the bottom plots 
give the results including the virtual ( terms. The solid curves are for 1PI kinematics, the 
dashed for PIM kinematics. 



the agreement between the 1PI and PIM kinematics further. This agreement also indicates that 
additional two-loop contributions not included in our expressions should be small. A similar 
trend is seen for the gg channel on the right-hand side of Fig. ^ The agreement between the 
NNLL 1PI and PIM scaling functions at low rj is significantly better than in the qq channel. Note 
however the significant divergence at large 77. Again, inclusion of the subleading contributions 
improves agreement over all 77. The improvement at larger 77, 77 > 0.1 is notable. We remark 
that the effect of the virtual ( terms is numerically small for both channels, in agreement with 
the arguments in Section IIIC of Ref. jHj. 

We now turn to our calculations of the hadronic total cross sections. We use the recent 
MRST2002 NNLO (approximate) parton densities [§| with an NNLO evaluation of a s . Our 
calculations use the exact LO and NLO cross sections with the soft NNNLL and virtual £ 
corrections and the full soft-plus-virtual scale-dependent terms at NNLO. In addition we mul- 
tiply the NNLO scaling functions by a damping factor, 1/^1 + 77, as in Ref. [Ij, to lessen the 
influence of the large 77 region where the threshold approximation does not hold so well. 

In Fig. |2J we present the NLO and approximate NNLO tt cross sections at VS = 1.8 TeV 
(left-hand side) and 1.96 TeV (right-hand side) as functions of top quark mass for fi = m. We 
also show the average of the two kinematics results. Both the 1PI and PIM cross sections are 
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Figure 2: The it total cross sections in pp collisions at y/S =1.8 TeV (left-hand side) and 1.96 
TeV (right-hand side) as functions of m for \i — m. The NLO (solid), and approximate NNLO 
1PI (dashed), PIM (dot-dashed) and average (dotted) results are shown. 



reduced due to the subleading terms. The average of the two kinematics results is just above 
the NLO cross sections for both energies. 

Our results using the CTEQ6M NLO parton densities ^U] are similar. At \/S = 1.8 TeV, 
averaging over the NNLO 1PI and PIM results with the two sets of parton distributions at 
H — m — 175 GeV, our best estimate for the cross section is 5.24 ± 0.31 pb where the quoted 
uncertainty is due to the kinematics choice. At y/S = 1.96 TeV our corresponding best estimate 
is 6. 77 ±0.42 pb. We note that differences with our previous estimates in jBJ[7j have as much to 
do with the use of the new parton densities as with the inclusion of the new subleading terms. 
We also note that the scale dependence of the NNLO cross section is negligible pQ. 
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